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Abstract. We propose an approach for quantum simulation of electron-phonon 
interactions using Rydberg states of cold atoms and ions. We show how systems of 
cold atoms and ions can be mapped onto electron-phonon systems of the Su-Schrieffer- 
Heeger type. We discuss how properties of the simulated Hamiltonian can be tuned 
and how to read physically relevant properties from the simulator. In particular, use of 
painted spot potentials offers a high level of tunability, enabling all physically relevant 
regimes of the electron-phonon Hamiltonian to be accessed. 

PACS numbers: 37.10.Jk, 37.10.Ty, 73.20.Mf, 32.80.Ee, 32.80.Qk 
1. Introduction 

Electron-phonon interactions lead to some dramatic effects in condensed matter systems, 
with lattice vibrations leading to phenomena such as superconductivity pQ and colossal 
magnetoresistance [2] . Even in the presence of overwhelming Coulomb repulsion, quasi- 
particle properties such as the effective mass can be strongly modified when electron- 
phonon interactions are strong. The importance of electron-phonon interactions has 
led to an industry in numerical simulation, but even very simplified models [31 H] are 
significantly more difficult to treat than other standard models in condensed matter 
such as Hubbard's model [5]. This difficulty lies in the very large Hilbert space that is 
needed to handle a potentially infinite number of phonons associated even with motion 
of a single electron. 

Recently, there has been a move to use systems of cold atoms as quantum simulators 
to investigate standard models of condensed matter physics [6j. The model that has 
been most successfully simulated is the Hubbard model of instantaneous local repulsion 
between electrons [S]. An optical lattice loaded from a Bose-Einstein condensate can be 
directly related to a Hubbard model, since there is a large energy penalty for two atoms 
to sit on the same lattice site [BJ. The use of cold atoms to simulate Hubbard models 
has led to direct observation of important phenomena such as the superfluid to Mott 
insulator transition [TJ [8]. However, there is no simple way of extending the approach 
to include phonons. 
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There are several classes of electron-phonon interaction that are of direct relevance 
to condensed matter systems. For polymers, which are very easy to deform, the 
hopping of electrons along the chain can be strongly modified by the presence of lattice 
vibrations. This leads to the Su-Schrieffer-Heeger (SSH) model, which describes a chain 
of CH groups in polyacetylene [3], which for a chain of N groups of mass M has the 
Hamiltonian, 

#ssh = - Yl ^+i>*( c I+i,<t ( V + H - c -) + y ~ + Y Yl ' ( X ) 

ia i i 

where the intersite hopping, t i+ i^ = to — g{ui + \ — Ui), c\ creates an electron at site i, g 
is the electron-phonon coupling strength, i is the hopping in the absence of vibrations, 
a is the electron spin, m is the group displacement and K is the spring constant for the 
phonons. The first term in this Hamiltonian is the electron kinetic energy from hopping, 
the second term is the potential energy from bond stretching and the third term the 
kinetic energy of the nuclei. In such strongly deformable materials, the vibrational 
amplitude of the atoms is large enough to affect the hopping integrals, since electrons 
typically hop more easily between atoms that are closer together. Thus the hopping, 
to, is augmented proportional to the relative displacement of neighbouring atoms. The 
energy spectrum for the second and third terms in the Hamiltonian, that relate to the 
phonon modes, may be solved in the absence of interaction by Fourier transforming, 
diagonalising and introducing creation and annihilation operators [U] . In this way, they 
are rewritten as ^2 ku f>^ku{dk v dkv + 1/2), where each phonon mode has energy fuxi kv . 
The combination d\ v d kv is the number operator for phonons of mode v with momentum 
k, where d hv and d ku respectively represent phonon creation and annihilation operators. 

General three-dimensional phonon fields are quantized in the usual way by 
substituting 

k,u 

with ^ ku the polarisation vector [9]. In the case of the standard SSH 
model, there is a single phonon mode, so the model is simplified and it, = 
Sfc \/h/2NMujk (dke~ tk Ri + d} k e tk ' Ri j . Creation and annihilation operators may be 

Fourier transformed using the relation, d\ = J2k d k e tk ' Ri / \fN , where d\ creates a 
phonon at site i. If the phonons are taken to be a dispersionless Einstein modes of 
frequency uo (which are a mean field approximation to the phonons, and are often a good 
approximation to optical phonon modes) then the momentum sum can be performed 
to obtain the expression, itj = y/H/2Mu Q [di + d\) . In the case of polymers, out of 
chain phonon modes may also be permissible, and if Einstein phonons are used as an 
approximation, the polarization vectors $, kl/ point in orthogonal directions. 

There are a number of other models of electron-phonon interactions in condensed 
matter. One of the most prominent is Holstein's molecular crystal model, which 
describes the type of interactions that can be found in many materials, where the local 
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electron density couples to local optical phonon modes which are usually associated 
with vibrations of other ions in the unit cell [3]. The Frohlich model describes the 
continuum limit of interactions between electrons and strongly polarisable materials, 
which originate when itinerant electrons generate a dipole moment in a medium [10J and 
can be generalised to lattice models [TT] (the latter are sometimes known as extended 
Holstein interactions). Holstein and extended Holstein interactions have the form, 

-^Holstein = — ^0 ^ ] C L C jo" 

ijo- j 

where g^ is related to the force between the electron and the ion, and for a standard 
Holstein model = gSij, where 6ij is the Kronecker delta. The local number operator 
for electrons of spin a at site i is ni a = c\ a Ci a , and djdj is the number operator for phonons 
at site j. The third term in the Hamiltonian represents a local Harmonic oscillator or 
Einstein mode. Since d^ + dj oc Uj where Uj is the atomic displacement, a simple physical 
picture arises: when an electron hops onto a site, an atomic displacement is formed, and 
as the electron leaves a site, the displacement relaxes. The continual displacement and 
relaxation of atoms generates a phonon cloud that moves with the electrons, increasing 
their effective mass. 

Even the simplest case of the electron-phonon model (the polaron) where a single 
electron moves on a chain has caused controversy since the polaron concept was 
introduced [TS]. In particular, the existence of self-trapping, where distortions of the 
lattice generated by the electron prevent the electron from hopping, has remained an 
open question. It is now known that self-trapping is not possible in the Holstein model, 
where the effective mass increases exponentially at strong coupling, but never becomes 
infinite [13J. On the other hand, there is tantilising evidence for self-trapping in the 
SSH model which may have consequences for polymer physics |14j . 

A number of schemes for quantum simulation of different interaction types 
between electrons and phonons have been suggested. A scheme for simulation of the 
Frenkel-Kontorova model to probe friction and energy transport has been proposed by 
Pruttivarasin et al. [15]. Bruderer et al. have proposed a system where an optical 
lattice is bathed in a BEC [TBI ELY] m which polaron effects have been observed [IB] . 
The use of cold polar molecules to obtain Holstein polaron effects has been discussed 
by Herrera et al. [IH] • It is also appropriate to note the use of Rydberg ions to simulate 
spin systems, where high-energy phonon states are used as part of the mapping onto a 
spin system. Porras and Cirac have made the mapping between a system of ions and an 
Ising model by making a Lang-Firsov transformation [20]. Their proposed system has 
the form of phonons coupled to a single Rydberg state via an interaction of the Holstein 
type, and the interaction is generated using a laser. Miiller et al. have proposed the use 
of trapped Rydberg ions to simulate spin chains [21] . Li and Lesanovsky have discussed 
the structural distortions associated with exciting a high-energy Rydberg state in a 
cold ion crystal [22]. We also note the proposed use of weakly dressed Rydberg states 
to simulate supersolids [23], strongly correlated gases [21] and collective many-body 
interactions [25J. 
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Figure 1. Schematic of the dressed Rydberg excitations. A set of cold atoms or ions 
is prepared. They are illuminated by two frequencies AE g s — and AE^p — Ap, 
where AE g s corresponds to the transition between the \g) and \nS) states, and AE^p 
to the transition between \h) and \nP) states. \g) and \h) are hyperhne levels and \nS) 
and \nP) Rydberg states with I = and I = 1 respectively, both with m; =0. The 
lasers are de-tuned from the transition so that the particles spend only a short period 
of time in the Rydberg states. This increases the lifetime of the particles by reducing 
the incidence of ionization, and also stabilises the particles against decay, although this 
reduces the effective dipole-dipole interactions between the atoms. Superpositions of 
the \h) and \nP) states move around the lattice as excitons. 

In this article, we present an approach to simulating the SSH polaron problem 
using Rydberg states of cold atoms and ions. We will show how a system of Rydberg 
atoms can be directly mapped to a Su-Schrieffer-Heeger electron-phonon Hamiltonian, 
which is of direct interest in the modelling of polymers and other condensed matter 
phenomena relating to highly deformable materials, and to our knowledge has not been 
investigated before in the context of cold atoms. The simulator is introduced and the 
mapping to corresponding models is carried out. We then discuss how to read properties 
of electron-phonon interactions from the simulator, and describe a simple experiment 
to investigate the self-trapping problem. 

2. Mapping between cold Rydberg atoms and a Su— Schrieffer— Heeger 
model 

We now consider the various terms in the Hamiltonian that describes the motion of 
excited Rydberg states, the interaction of the Rydberg state with ions and atoms, 
and the interaction of atoms and ions with each other. We begin by discussing how 
electron-phonon interactions can be simulated in a system of cold Rydberg atoms such 
as rubidium in a deep optical lattice where there is a single atom per lattice site. The 
deep optical lattice is essential in this case so that hoppings of atoms between sites are 
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suppressed and there is no double occupancy of atoms. Systems of cold ions will be 
discussed later. 

Following the prescription in Ref. [26], it is possible to set up a stable (long lifetime) 
Rydberg system where the excitations are superpositions of low-energy hyperfme and 
high-energy Rydberg states. Such a system is initialized with 2 hyperfme states \g) and 
\h), so that (for example) on a chain the system is in a state — \ ■ • ■ ggghggg ■ ■ ■). 
Then two lasers detuned from the \g) -H- \nS) and \h) <H- \nP) transitions of the 
cold atoms from the ground to Rydberg states \nP) and \nS) are switched on so 
that all but one atom is in the dressed state |0) = cs\nS) + c g \g). Excitons in state 
|1) = cp\nP) + Ch\h) hop in the system, as illustrated in figure [TJ Here cs,c g ,cp and 
Ch are constants. As described in Ref. [26] the van Vleck perturbation theory (which 
acts to project an interacting Hamiltonian onto an effective hopping Hamiltonian) can 
be used to obtain an effective Hamiltonian for the motion of excitons, 

H eS = J2(E 2 + E, + W) |^>(^| + ( 4 ) 

where the wavefunction \ jti) represents a single exciton at site i. There are onsite energy 
offsets as a direct result of the dressing of the Rydberg states, which are calculated to 
2nd order as, 

E 2 = (N — l)a 2 s A s + apAp (5) 
and to 4th order in the van Vleck scheme as, 

E 4 = -{N - l)a 4 s A s + a P A P - (N - l)a 2 s a 2 P (A s + A P )} (6) 
and there is an additional energy shift due to dipole-dipole interactions of, 



W = a 2 s a 2 P 



A hopping term between exciton states on sites i and j also emerges from the 
perturbation theory, 

% = ^Wpt%2 (8) 

where Vij = Vij/{Ag + A P ). N is the number of ions in the chain, as = VLs/2As 
and ap = Qp/2Ap are dressing parameters. and Ap are the detunings of the 
illuminating radiation away from the frequencies corresponding to the g-nS gap, AE g $, 
and the h-nP gap, AEhp, respectively. The Rabi frequencies are denoted fls and Qp. 
The undressed interaction between Rydberg states on different ions is given by the 
dipole-dipole interaction, 

where \x is the effective dipole moment on the Rydberg atom and Rk is the vector to 
the kih atom in the optical lattice. 
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All states are set up with mi = 0, and the atoms are illuminated with light which is 
polarized out of the plane, fixing the dipole moments perpendicular to the interatomic 
vectors. With this choice of laser polarisation, transitions to the \nP),mi = 1, —1 states 
are suppressed |27J. A further utility of this set-up will become apparent when cold 
ions are discussed. It is worth noting that the effective Hamiltonian can also be written 
using the second quantised notation as, 

H eS = S £ J {E 2 + E i + W) a\ ai + % a \<*j (10) 

where a\ creates an exciton and annihilates an exciton at site i. In the following, the 
excitons will take the place of the electrons in the SSH and Holstein models. Note that 
the excitons are hardcore as there cannot be more than one exciton per site (there is 
either an exciton associated with a site, or there is not). 

Terms in the Hamiltonian Q have the following origin: The first term (E2 + E4) is 
simply an onsite, constant offset of the energy which would be absorbed into a chemical 
potential. It arises because of the interaction between the lasers and the atoms, and 
would be present even in the absence of dipole-dipole interactions between the atoms. 
The energy offset term W arises because the local dipole moment induces moments in 
other atoms, which then cause an energy shift by re-interacting with the original dipole. 
The most important term here is the long range hopping term Vij which is the origin of a 
Su-Schrieffer-Heeger interaction as discussed below. The corrections to this interaction 
also have their origin in induced dipoles. 

The following discussion is for arbitrary dimension and works equally well for optical 
lattices in planes and on chains, and for phonons that can vibrate on a line, plane or 
within a full 3D space. Phonons can be introduced to the system by displacing the 
atoms from their equilibrium positions. Atomic displacements within the optical lattice 
are independent as no interaction is carried through the optical lattice. Therefore, the 
vibrations of the atoms in the lattice are a faithful representation of Einstein phonons. 
One interesting possibility for 2D or 3D vibrations is to distort the harmonic potential 
of the optical lattice to lift the degeneracy of the phonon modes. The Hamiltonian for 
these modes will be, 

H ph = J2 huov(dt k d vk + 1/2) = ^M vi d vi + 1/2). (11) 

uk ui 

The introduction of small in-plane phonon displacements, itj, at site i causes the 
interaction between Rydberg states to become, 

Vkl = -fs-r — ^ is' ( 12 ) 

\R k + u k -Ri- ui\ 6 

which can be Taylor expanded if the displacement is small, noting that it is only 
necessary to expand to 1st order in u as this is the lowest order for non-trivial interaction 
with phonons. This leads to: 

V ( R+U ^\w~^m^- (13) 
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This expansion can be applied to the off-site terms in equation [8] to obtain, 



V tJ (u) 



a 2 s a 2 P jj 2 



1 - Ji A /R 



VI 



QtLij ' R j ij 



R 3 - 



3itjj ■ Rij 



I? 



R 10 



+ 0{u 2 



(14) 



where /z 2 = /i 2 / (A s + A P ) and ity = Ui — Uj. 



Quantisation of the atomic displacements using equation [2] leads to an extended Su- 
Schrieffer-Heeger (SSH) model which is general for both cold ion or cold atom systems, 



H, 



ext-SSH 



a s a P jj, 



n A /R% Rl 

dkA e 



ik-Rj 



1 - 



h 



R^ ^ V 2NMu ku 



Ri 



ij-^ku 

e -*** Jlt ] +<$ hv [e ih ' R i -e tk - R ']) 



(15) 



Here, af/a$ create/annihilate excitons on site i and M is the mass of the atoms (or ions). 

A further simplification can be made for atoms in the optical lattice because the 
oscillators are independent and local, thus there is no momentum dependence in the 
phonon modes and Uk u = Wqv and ^ ku = £ 0u so that the momentum sum only applies 
to the phonon creation and annihilation operators. Moreover, the set of £ 0u is now 
orthogonal. Fourier transforming these operators, the local nature of the model with 
Einstein modes becomes clear, 



H. 



ext-SSH 



E+ Ctc 



a 2 s a 2 p/j, 2 



tf/R% R% 



(16) 



l 



This Hamiltonian is an extended version of the Su-Schrieffer-Heeger model. Note that 
there can still be several phonon modes since even if the atoms are arranged in a chain 
or sheet, the harmonic traps can be three dimensional. Making the comparison with the 
SSH Hamiltonian (Eq. [TJ, it is possible to assign values to to and g. The form of the 
coefficients mean that the relative strength of the parameters can be tuned by changing 
the inter-ion spacing in the lattice since the Hopping term, t, has 1/R 3 functionality, 
whereas the electron-phonon coupling term g ~ 1/R A . The phonon frequencies can be 
changed by altering the depth of the optical lattice, giving full control over all parameters 
in the model. We will discuss this tuning in more detail in section |4| 

Making a similar Taylor expansion, phonon interactions via the on-site term in 
equation [7] can also be considered, 

1 1 

(17) 



E 

k^j 



-^/R% + Qu kr R k f^/Rl 3 



Er 

k+i 



ft'/Rt 



6u k j ■ R k jH 

RU^W/Rt 



fcjj 
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interaction 



Figure 2. System of cold Rydberg ions / atoms used to simulate a polaron, annotated 
with Hamiltonian terms. 



so that the local Hamiltonian becomes, 

#ext-Hoi(«) = ^ a\ai{A s + Ap)a 2 s a 2 



p 



A ^ 4 /^L 



(19) 



the first term in this expression leads to a constant energy term. Once quantized, the 
second term leads to an extended Holstein type electron-phonon interaction. However, 
the prefactor 1/R 7 indicates that this term will be very small in comparison to the SSH 
style term. 

Combining all Hamiltonian terms from the preceding discussion, the following 
Hamiltonian is found (up to a constant term): 

H = -ffext-SSH + -Hext-Hol + -^ph (20) 

A schematic of the simulator with the effective interactions is shown in figure [2] 

In order to gain insight into any differences that long range hopping make to the 
extended SSH Hamiltonian, we make perturbation theory calculations for three cases: 
The standard SSH model, an extended SSH model with 1/R 3 hopping and the extended 
SSH model with both 1/R 3 hopping and 1/R 4 interaction strength. In all cases, the 
model is on the chain. A single Einstein phonon mode in one-dimension is considered. 
The lowest order perturbative contribution to the self-energy of the polaron at absolute 
zero is: 

A f W , 9k,k+q 



E(w,fe) = — / dq (21) 

2tt J- n u -u - e k+g + id 

where an integration over the energy has been performed for the phonon propagator. 
The Fourier transform of the interaction vertex is, 

OO j 

g Kk+q = 2ig V [sin(n(fc + q)) - sin(nA;)] (22) 
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Figure 3. Effect of long range hopping and vertex on the polaron shift of the SSH 
polaron formed from interaction with one Einstein mode. Here A = g 2 /woio is the 
dimensionless electron-phonon coupling and ujq = O.li. There is no qualitative change 
in the polaron shift as long range effects are switched on. 




Figure 4. Effect of long range hopping and vertex on the effective mass of the SSH 
polaron formed from one Einstein mode. 



and the non-interacting dispersion, 

E x cos(nk) , . 

„, <w (23) 

The dimensionless electron-phonon coupling A = g 2 /uoto and 5 is a small real number. 
From this, the polaron energy can be calculated through the expression, 

w-e fe -Re£(u;,A;) = (24) 

and the effective mass from the second derivative of the polaron energy with respect to 
momentum, k. 

Results from the perturbation theory are shown in Figs. [3] and |4j The polaron 
shift shown in Fig. [3] is the difference in energy between the interacting and non- 
interacting cases. For all three cases, the form of the polaron shift is qualitatively 
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similar. Introduction of long range hopping decreases the polaron shift as the polaron is 
essentially less tightly bound to a site and hopping becomes easier, whereas introducing 
the long-range interaction moderately increases the shift because the effective interaction 
is slightly stronger. A stronger change can be seen in the effective mass as the hopping 
range is increased (Fig. [4]) but the core physics remains the same. The divergence in 
the effective mass signifies the breakdown in the perturbation theory. 

We note that the Rydberg excitation is Bosonic. For polaron problems, which have 
a single electron, this is no problem because the electron has no other particle with 
which to exchange. It is also possible to simulate the ground state of a bipolaron (two 
interacting polarons) since the ground state of two interacting Fermions is always a 
singlet with symmetric spatial component of the wavefunction, with the antisymmetry 
of the wavefunction imposed by the spin degrees of freedom j2H] ■ Thus, most observables 
in the absence of a magnetic field (e.g. total energy, effective mass) are identical in the 
two Fermion or two Boson ground states. Since there can only be one exciton per site, 
the Bosons are hard core. 



3. Cold Rydberg Ions: Interactions with acoustic phonons 

An alternate possibility for tuning the Hamiltonian is to use cold ions instead of cold 
atoms. For a system of cold ions confined to ID, the ions form a chain, held in place 
by the direct Coulomb repulsion between the ions. Cold ions confined to 2D will form a 
Wigner crystal with a triangular lattice. The phonon modes have been calculated in the 
harmonic approximation for the triangular lattice [23 ED] and have an acoustic nature. 



Note that these modes are not independent and the full Hamiltonian 15 applies here. 
The ions considered here are sufficiently heavy to justify the harmonic approximation. 
Therefore, by changing from rubidium atoms to strontium ions, systems with acoustic 
phonons are accessible. 

Interaction between the Rydberg states and phonons occurs where there are 
interactions terms between the dipole moments. For a system of cold Rydberg ions 
with small dipole moment, the interaction potential is 



, . „ , 1 (IL — Ra) ■ (n — r,-) 
V{R U R 3 , n, rj) = —+ y — J -j-^ 3 ± 

rj - 3(n tJ ■ Tjf + rj - 3(n tJ • r 3 f 
+ 2^ 

The system should be prepared in the state discussed at the beginning of section [2j 
which simplifies the interaction to, 

_ . 1 u 2 
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Figure 5. Relative strength of on-site and near-neighbour energies in an optical 
lattice. f2 = 10MHz and A = 100kHz. The polarons formed from weak coupling with 
high-energy phonons can be investigated. 



The monopole ion-ion interaction is only responsible for the phonon modes and the 
effective dipole-dipole interaction is identical to the cold atoms case. Thus, beyond 
a change in the form of the phonons, the resultant exciton-phonon interaction is the 
same for atoms and ions. We note that a fully excited \nS) ion moves in a different 
potential to an ion in a fully excited \nP) > state. Since the ions are dressed, they have 
only a very small admixture of the Rydberg states, so the motion of the ions is mainly 
determined by the ground state properties of the ions which are very similar for \g) and 
\h) states, and therefore it is safe to assume that all ions move in approximately the 
same potential. 

4. Experimental considerations 

So far we have discussed the Hamiltonian of the atomic/ionic system and its condensed 
matter analogues. Now we sketch the manner in which experiments may be performed. 

In the case of atoms or ions, we assume that each lattice site is occupied by a single 
atom/ion, and each site may be addressed individually in such a way as to allow the 
preparation of the initial simulator state \g) or \h) on each atom/ion. We assume also 
that the state of each atom may be measured by use of a laser tuned close to resonance 
with a closed transition involving the ground state of interest. Thus, the state of each 
atom/ion is detected in fluorescence. 

In the case of atoms, the system may be prepared from a Mott insulator phase 
where a single atom occupies each lattice point. Alternatively, there exist proposals 
to obtain single occupancy in non-degenerate (and hence much lower density) atomic 
ensembles [31J . For ions at low temperatures a suitable lattice is necessarily formed; it 
is possible to cool the motional states of the ion crystal using finely tuned lasers |32j . 

First, we discuss the Rydberg states which cause the atom- atom interaction in our 
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SSH model. We couple our two hyperfine ground states \g), \h) to the \nS) and \nP) 
Rydberg state, respectively. In the following we discuss the atom Rb 87 . These particular 
Rydberg states interact according to the dipole-dipole Hamiltonian in Eq. [9| In detail, 
this takes the form: 

V= (27) 

where the vector connecting the two atoms is R, d n s- n p the dipole matrix element 
between the states \nS) and \nP) and Y2 tm - m >(R) is a spherical harmonic. The lasers 
select only the m = and ml — magnetic sublevels. In general, the dipole-dipole 
interaction scales with n A / R 3 and the lifetime of the Rydberg states scales with n 3 , which 
together mean that choosing the highest energy Rydberg state possible leads to the best 
possible properties of the atoms and ions for the simulator. However, experimentally we 
stress that for a given laser intensity, the Rabi frequency of the dressing laser scales with 
7?r 3 / 2 , and the shift of the \nS) atoms due to their mutual interaction (the "blockade 
shift") scales with n 11 / R 6 . Thus one cannot go to arbitrarily high n without shifting 
the \nS) states too far from the dressing laser, or by making unreasonable demands 
on the laser intensity. We choose the n = 48 states for which the laser power required 
to drive the \g) — > \nS) and \h) — > \nP) transition should be reasonable (the Rabi 
frequency for these transitions scales with n~ 3 ^ 2 ), but where the 48S - 48P interaction 
is 1252 MHz//im -3 ; large enough for our purposes. The lifetimes of these states are: 
7"48S = 106/is and r 48 p = 117/is, but the effective lifetime of the dressed states are 
lengthened by a factor a^ 2 and a J, 2 . Using this expression for V and reasonable values 
for the Rabi frequencies Q$ — &>p = 10 MHz and detunings As = Ap=100 MHz we 
calculate the nearest neighbour exciton hopping rate for a separation of 2.5 micron 
to be 0.19 KHz, and a dressed state lifetime of 48 ms ~ 54 exciton hopping times. We 
note that for these choices of Ag,Ap and fls, &p, diverges below a separation of 
about 2.2/iin, which sets an additional constraint on the types of experiments we can 
perform. We also know that the \nS) atoms will interact with each other due to the 
van der Waals interaction V v dw = Cq/R^, where Cq is the appropriate van der Waals 
coefficient. For our 48S* states, C 6 « 1050 MHz//mi 6 . 

Now we turn to considering some of the experimental values required for modelling 
an SSH Hamiltonian. We have in mind two set ups which illustrate the utility of the 
model and some of the experimental issues that must be addressed. In the first case, a 
linear chain of atoms (of mass M) is confined by a far-detuned optical lattice formed from 
a pair of intersecting laser beams with wavelength A. As mentioned above we stipulate 
that there be a single atom per lattice site. With counter-propagating lasers, such a 
system is characterised by the recoil energy E r = 2iih 2 /MX and the potential depth V , 
which depends on the laser intensity. Assuming that the atoms are trapped near the 
potential minima, the phonon frequency will be Uo = 2E r ^Vo/ E r /h, the atomic hopping 

energy J = (4 / y/ir) E r ( Vo / E r ) 3 / A e^ Vo ^ Er and we estimate the on site atomic repulsion 
energy as U = ^8/7r(a s /a )E r (Vo/E r ) 3 / 4 } with a s the s-wave scattering length for the 
atoms - 5.3 nm for Rb 87 , and Qq the length of the ground state wavefunction. In order to 
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initialise the system, it is necessary that the atoms be individually addressable, which 
determines the minimum lattice wavelength to be around 1 /im, however we adopt the 
more conservative 2.25 /im lattice wavelength considered above. In order that atomic 
tunneling rate J be suppressed, we require that Vq » 20E r = 4.53 kHz; for the 2.25 
jum lattice spacing, Vy = 1.43 kHz and Uq = 2.03 kHz. We then have J ~ 0.001 kHz, 
U = 2.3 kHz, increasing the localisation further. We note that this phonon frequency 
is a large fraction of the site depth Vq. Also, for atomic spacings of 2.5 microns and 
upwards, V v dw < 24 MHz, smaller than, but comparable to, the detunings As and Ap 
we have chosen for our example. 

The tunability of the cold atom system in an optical lattice where the atomic 
potential Vi att = Vq sin 2 (/cx), Vq = 20E r is illustrated in Fig. [5] as a function of the 
lattice wavelength. The figure shows the relative energies associated with the trap 
depth, the phonons, the nearest-neighbour exciton-phonon coupling and the exciton 
hopping. Using realistic parameters, the phonon energy and coupling energy can be 
swapped. The resulting SSH model is in the weak coupling, high phonon frequency 
regime. 

In the next experimental scenario we consider how a far more tunable system can 
be constructed by using painted potentials [33]. In this system the lattice potential is 
an array of Gaussian laser spots, with power P and focused to a waist Wq, in which case 
the traps have a potential V irav = Vb pt e~ 2r2//u, o where the trap depth is Vopt = ZP/kWq 
and the phonon frequency is u = (1/itfo) v Vopt/M. The trap spacing, d, may be chosen 
such that tunneling is negligible by choosing d > wq, and the phonon frequencies may 
then be determined by the laser power P. In this case, we need to choose higher Rabi 
frequencies and detunings Qs,p — 20 MHz and A$,p = 100 MHz respectively. For 
these parameters, and with the larger atomic separations (> 4/xm), we calculate that 
the dressed state lifetimes are ~ 12 ms, which is about 85 Th op , and that the blockade 
shift of the nS states is less than 2MHz. We plot the the parameters associated with 
this potential in Fig. [6j Here, all hierarchies of phonon energy, coupling constant and 
exciton hopping can be achieved. A major advantage of the painted potentials is that 
much deeper potentials can be achieved, but phonon frequencies can be much smaller, 
putting the effective model in the physical regime of a condensed matter system. Also, 
rings of atoms can be created to explore periodic boundary conditions. 

For modelling the SSH model using cold ions, we have not undertaken any detailed 
calculations to assess our simulator proposal, but we can make some points here. Ions 
trapped in an ion trap have phonon frequencies in the 100 kHz - 2 MHz range, and 
spacing on the order of several /im. Thus, the phonons are very much higher in energy 
than the atomic case. Recent studies of Rydberg states of ions has shown that the 
interactions between the ions are very weak - additional complications involving RF 
dressing of the rydberg atoms are required [34] . This means that implementation of 
the cold ion extension to simulate polarons formed from acoustic phonons in regimes 
relevant to condensed matter physics will be a major experimental challenge. 

We end by describing an experiment to explore the self trapping regime of the SSH 
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Figure 6. Relative strength of on-site and nearest-neighbour energies in painted 
potentials. All hierarchies of energies can be realised. The parameters here are 
detuning of A = 100kHz, Rabi frequency f2 = 25MHz and the waist is 2.5/im. (top) 
The painted potential is made with a P = 0.75mW laser, and the exciton hopping and 
coupling energies cross above the phonon energy so that the regimes where phonon 
energy < hopping < coupling and phonon energy < coupling < hopping can be 
explored, (bottom) P = 0.75mW, and the exciton hopping and coupling energies 
cross above the phonon energy so that the regimes of phonon energy > hopping > 
coupling and phonon energy > coupling > hopping can be explored. 

model which is a subject of current interest |14j . We consider studying the phenomena 
of self-trapping in the SSH model on a circular array of Gaussian spots - that is, a ring 
of traps. A particular lattice point is chosen for the location of a single excitation, and 
then the dressing fields are rapidly switched on. Since this impulse contains multiple 
momenta, the localized exciton state will expand into the lattice. The result of the 
simulation may be obtained by extracting the location of the excitation as a function of 
time. This may be achieved by rapidly switching off the dressing potentials, and then 
measuring the state of the atoms/ions at each lattice point. In the self-trapped regime, 
the exciton hopping will be suppressed which will be observed as a narrowed distribution 
of final exciton locations. The hopping leads to a quantum random walk [35J so the mass 
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can be extracted from the spatial distribution of an ensemble of measurements. Note 
that there are significant differences between the distributions arising from classical and 
quantum random walks. 

5. Summary and outlook 

In this article, we have shown how systems of cold Rydberg atoms and ions can be used 
as a simulator for electron-phonon interactions, and derived the mapping onto a Su- 
Schrieffer-Heeger model, which is used to simulate strongly deformable materials. We 
have also discussed how readout of physically relevant properties could be carried out. 
Further, we have discussed the parameters of the model that can be achieved, showing 
that rubidium atoms in painted potentials can be used to access physically relevant 
couplings, phonon frequencies and exciton hoppings. 

It is relevant to discuss scalability here. A drawback of the current approach is that 
the excitations are Bosons rather than Fermions. It may be possible to simulate full 
fermion-phonon interactions using a more complicated two layer 2D optical lattice with 
dressed Rydberg atoms in the lower layer and an incomplete layer of Rydberg atoms in 
the upper layer. The Rydberg atoms in the upper layer will interact with those in the 
lower layer through the dipole-dipole interaction, leading to a Hubbard-Frohlich model, 
but with the caveat that the system is much more complicated [US]- We note that much 
has been learnt by simulating Bose-Hubbard models with cold atoms, so we consider 
that the lack of Fermionic excitations is currently a minor drawback. We hope that 
the proposed system will stimulate experimental work in the area, and encourage the 
development of quantum simulators for electron-phonon interactions which work with 
Fermionic excitations. In the case of polymer physics, the quantum simulations are 
relevant as they may shed light on ongoing debates about self-trapping phenomena [Hj . 
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